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ABSTRACT 

Aims. Very high energy (VHE) gamma-ray emission from a distant source (z > 0.2) can be efficiently absorbed my means of the 
electron-positron pair creation process. Analyses of the unabsorbed spectra imply that the intrinsic TeV emission of some blazars is 
hard, with spectral indices 0.5 < a < 1. The absorption depends on the level of extragalactic background light (EBL) that is difficult 
to measure directly. This implies that it is difficult to estimate the slope of the intrinsic TeV emission. To test our blazar emission 
scenario that is capable to reproducing the hard spectra, we therefore used the observations made by the Fermi Gamma-ray Space 
Telescope in the unabsorbed MeV-GeV energy range. 

Methods. We assume that the X-ray and gamma-ray emission of TeV blazars is produced in a compact region of a jet uniformly filled 
by particles of relatively high energy (y > 10 3 ,£ = ym c c 2 ). In other words, we assume a low energy cut-off in the particle energy 
distribution. The emission produced by the particles with this energy spectrum can explain hard intrinsic spectra in the energy range 
from MeV up to TeV. We demonstrate how to estimate the basic physical parameters of a source in this case and how to explain the 
observed spectra by a precise simulation of the particle energy evolution. 

Results. To test our estimation methods, we use the observations of two blazars with exceptionally hard spectral indices (a < 0.5) 
in the MeV-GeV range and known redshifts: RGB J0710+591 and 1ES 0502+675. The estimated values of the Doppler factor and 
magnetic field are compared with our numerical simulations, which confirm that the particle energy distribution with a low energy cut- 
off can explain the observed hard spectra well. In addition, we demonstrate that the radiative cooling caused by the inverse-Compton 
emission in the Klein-Nishina regime may help us to explain the hard spectra. 
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c-3 1. Introduction 



The emission of some blazars is observed from radio frequen- 
cies up to VHE gamma rays. Collating observations from differ- 
ent energy ranges, we can show that the spectra of these objects 
contain two characteristic peaks (in vF(v) plots). The first peak 
appears in the X-rays from a few keVs to a few hundred keVs 
and the second peak is observed around an energy of a few TeVs 
(e. g. Ghi sellini, G. et al. I19981 Massaro et al. l2004l N ieppola et 
al. 120061 Ghisellini, G. et al. 120101 Abdo et al. l2010bt . Since the 
TeV emission is the most prominent feature in those sources we 
often call them TeV blazars, which constitutes a relatively small 
group of sources within the blazar family. The gamma-ray emis- 
sion of most blazars peaks in the MeV-GeV range. Moreover, 
TeV gamma rays are efficiently absorbed by the EBL. This in 
addition limits the number of observed TeV blazars. 

The high energy emission of blazars is believed to originate 
inside their jets. However, the observed variability on time scales 
from days (e.g. Catanese et al. 1 19971 Fossati et al. 2008 ) down to 
a few minutes (e.g. Aharonian et al. [2007a) indicates that only 
a small part of a jet is radiating at high energies. This kind of 
emission requires high energy particles, which can gain energy 
by acceleration at the front of a shock wave inside the jet. This 
is the simplest explanation of the acceleration process, which 
assumes that a small fraction of the jet bulk kinetic energy is 
transferred to the particles. In the presence of a magnetic field, 
the particles can radiate this energy through synchrotron emis- 
sion, generating the first peak in the spectrum. Some fraction of 



the synchrotron photons can then be up-scattered to higher ener- 
gies. This is the inverse-Compton (IC) scattering that gives the 
second peak in the TeV energies. This simple scenario has been 
proposed many times as part of models of the VHE emission of 
blazars (e.g. Dermeret al. 119971 Bloom & Marscher 1996, Inoue 
& Tak ahara |1996l Mastichiadis & Kirk |19971 Katarzyriski et al. 
1200Tb . 

A TeV gamma ray photon travelling through the intergalac- 
tic medium can interact with an infra-red photon producing an 
electron-positron pair. This causes absorption of the emission 
above a few hundred GeVs. To calculate this absorption, we have 
to determine the level of the intergalactic radiation field, which 
is difficult to measure directly but can be estimated from sim- 
ulations of star light production in evolving galaxies. There are 
many different solutions for the level of the absorption (e.g. De 
Jager & Stecker 120021 Keneiste et al. 120041 Franceschini et al. 
2008 Kneiske, & Dole l2TjT0l l. Thus the slope of the intrinsic TeV 
emission cannot be accurately calculated. Nevertheless, the TeV 
emission of some relatively distant sources (e.g. 1ES 1101-232 
z = 0.186, Ahar onian e t al. 120061 or 1ES 0229+200, z = 0.14, 
Aharonian et al. 2007b ) indicates that the intrinsic spectra are 



hard {F 



with 0.5 < a < 1), even if we assume as low as 
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possible absorption. 

It is difficult to explain how these hard spectra could be cre- 
ated. Most of the standard emission models assume a power-law 
or broken power-law particle energy distribution with a slope 
n 2 (where the number of particles N oc y~ n ) and a mini- 
mum energy equivalent to y m 1 . This gives the spectral index 
of the synchrotron emission a = (n - l)/2 = 0.5. This emis- 
sion is up-scattered to higher energies by the same population of 
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the particles, hence the index of the inverse-Compton spectrum 
should be similar to the index of the synchrotron emission. Note 
that high energy particles (y > 10 5 ) do not scatter efficiently be- 
cause of the Klein-Nishina (KN) restrictions. The broken power- 
law particle distribution (n\ and «2 below and above the break, 
respectively) implies that there is a broken synchrotron spectrum 
(a\ and o^), for which the spectral index in the KN regime is ap- 
proximately a =a 2»2 — a\ = 2.5 for «2 = 4 (Tavecchio et al 
[19981 

A question to ask is how we can explain the hard spectra 
from MeVs up to TeV gamma rays. One of the simplest solu- 
tion is to assume a cut-off in the low energy part of the electron 
spectrum, which basically means y m » 1. This solution was 
proposed for the first time by Katarzynski et al. (2006a) to ex- 
plain the observations of 1ES 1101-232 in the TeV range. This 
cut-off leads to an additional break in the synchrotron spectrum 
between IR and X-ray energies. The spectral index of the first 
part of the spectrum, from radio up to X-ray energies, becomes a 
constant ao = -1/3. This is a "tail" of the synchrotron emission 
produced by the low energy particles. This part of the radiation 
field with such a hard spectral index can also be up-scattered. 
Therefore, the IC spectrum can also be very hard up to limiting 
value a = -1/3. Simultaneous optical and X-ray observations of 
1ES 0229 +200 performed by the Swift satellite (Tavecchio et al. 
2009 ) shows an abrupt break between the optical and the X-ray 
range. This strongly implies that there is a low energy cut-off to 
the particle energy distribution. 

In the present work, we explore this idea of a low energy 
cut-off. We focus on the unabsorbed MeV-GeV range where the 
spectral index of the intrinsic emission is observed directly. We 
demonstrate how to estimate the basic parameters of a source 
and how to explain the observed spectra by a simulation of the 
particle energy evolution. We apply our estimations and mod- 
elling to the X-ray and gamma-ray observations of two blazars 
with known redshifts. 



2. Observations 

The spectra of TeV blazars in the MeV-GeV range are always 
hard with spectral indices a < 1 . The average expected value of 
the spectral index in this particular range is a ~ 0.5. In this work 
we focus on the sources with a < 0.5 that we call exceptionally 
hard. 

For all blazars observed by the Fermi Gamma-ray Space 
Telescope (Atwood et al. 120091 from 20 MeV up to 300 GeV, 
only in some cases were spectra published in the Fermi Large 
Area Telescope First Source Catalog (Abdo et al. 120 1 0al > excep- 
tionally hard with the photon index T < 1 .5 (where the photon 
index is greater by unity in comparison to the spectral index). 
Moreover, in many cases we either do not know the distance to 
the observed sources (e.g. RBS 158 r = 1.34 ± 0.33, RBS 621 
r = 1 .42 + 0.32) or the spectra were obtained with a relatively 
large uncertainty (e.g. 1ES 1 101-232 T = 1.36+0.58). Therefore, 
we selected only two objects for our test. 

The observational parameters of the first object RGB 
J0710+591 (z=0.125) in the Fermi catalogue are Y = 1.28+0.22, 
and F = (1.46 ± 0.49) 10" 11 erg crrT 2 s _1 , where F is the en- 
ergy flux in the 100 MeV to 100 GeV range. This source was 
also discovered by the VERITAS gamma-ray observatory in the 
TeV range (Acciari et al. 1201 Oi l. The photon index obtained by 
VERITAS above ~ 300 GeV was F = 2.69 + 46. In our analysis, 
we also use X-ray observations obtained by the Swift satellite 
in the energy range from 0.2 to 10 keV (Acciari et al. 1201 01) . A 



notable characteristic of this object is that the X-ray peak is ob- 
served at the level a few times higher than the GeV- TeV peak. 
This indicates that the synchrotron cooling is dominating in this 
source. 

The photon index of the second source 1ES 0502+675 
(z=0.341) in the Fermi catalogue is rather soft F = 1.75 ± 0.1 1. 
However, the MeV-GeV spectrum of this source is quite com- 
plex. Owing the limited photon statistics, the observations pub- 
lished by Abdo et al. (2010E) and Abdo et al. (1201 Pel ) seem 
to display a broken power-law spectrum of IC emission. This 
is rather unexpected and not predicted by any model of blazar 
emission. However, this may have a simple explanation if we 
assume for example that two independent jet components were 
observed at the same time. The break in the spectrum appears 
around an energy of 1 GeV. Fitting a single power-law func- 
tion to the observations above this energy, we obtained F = 
1.35 + 0.20 and F = 1.46 + 0.18 for the data sets above quoted, 
respectively. The VERITAS team reported the detection of this 
source in the TeV energy range (Ong 2009). However, no public 
information about the observed flux and slope has been publicly 
available so far. In addition we use quasi-simultaneous observa- 
tions from Swift (Abdo et al. 2010b). The comparison between 
the X-ray peak level and the GeV observations shows that in this 
particular source the radiative cooling generated by the IC emis- 
sion is the predominant cooling process. 



3. Estimations 

To estimate the basic parameters of a source of VHE emission, 
we adapted approach proposed by Tavecchio et al. (1998). We 
assume a spherical source of radius R, filled uniformly by a con- 
stant magnetic field (B) and relativistic particles {K - density). 
In addition, we assume that the particle energy distribution is a 
broken power-law (ti\, n%, y\, - break energy) and that the source 
travels with relativistic velocity (6 - Doppler factor) inside a jet. 
Finally, we included a low energy cut-off in the particle energy 
spectrum (y m - minimum energy). This is our main extension 
to the initial approach that let us estimate the basic physical pa- 
rameters from the exceptionally hard spectrum. Even this sim- 
ple scenario has eight free parameters. Therefore, using obser- 
vational quantities we can only estimate the relation between the 
main parameters B and 6. However, we can estimate this relation 
in four different ways. 



3.1. Peak positions 

The simplest estimation of the B - 6 relation involves the com- 
parison of the synchrotron peak position (v s ) with the inverse- 
Compton peak frequency (v c ). The positions of the peaks do not 
depend on y m thus we can directly use the formula derived by 
Tavecchio et al. JT9981 



B = (l+z) 



2.8 x 10 6 v c 



(1) 



In Fig.[T| we illustrate how an approximate IC spectrum evolve 
with time for different values of y m , from y m > 1 up to jt,. The 
spectra are calculated according to the prescription given in the 
appendix. This simple test shows that the IC peak remains at the 
same position for the wide range of y m values with the exception 
for y m =i jb and this gives the limit to the above formula. The IC 
spectrum with a - -1/3 that could be obtained for y m - jb was 
never observed directly. The observed spectral indices described 
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approximated synchrotron spectra 
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Fig. 1. Approximate synchrotron and IC spectra calculated for 
different values of y m , from almost unity up to jf,. This test 
demonstrates how the spectral index in the MeV-GeV range de- 
pends on the y m value. The dashed line shows the difference be- 
tween the estimated (a\ = 0.5) and calculated {a 0.65) spec- 
tral slope below the peak. 



in the previous section are in range from 0.3 to 0.5. Therefore, 
we can safely use this B and 6 relationship for our estimations. 

The spectra in Fig. Q] were calculated for n\ — 2 and «2 = 4, 
which directly gives the spectral indices of the synchrotron emis- 
sion below (ari = 0.5) and above (aj - 3/2) the peak, respec- 
tively. The value of n\ = 2 is rather typical for many differ- 
ent astrophysical sources. This may be the result of the first- 
order Fermi acceleration at a shock front or the synchrotron or 
the inverse-Compton cooling in the Thompson regime. This in- 
dex is crucial for the emission of TeV blazars because this part 
of the particle spectrum provides synchrotron photons for the 
IC scattering. The IC spectrum in the MeV-GeV range comes 
mostly from the scattering in the Thompson regime. Therefore, 
the spectral index of this emission should be similar to the index 
of the scattered radiation field. However, the complexity of the 
scattering (see Fig. lA.U produces a curved spectrum below the 
peak. In the MeV-GeV range, just below the peak, the spectrum 
can be approximated well by a power-law function with the in- 
dex a =i 0.65. This differs significantly from the estimated value 
a = 0.5 that is clearly illustrated in Fig. Q] Note that this con- 
cerns only the case where y m ~ 1 because y m controls the slope 
in the other cases. Since the somewhat classical value of n\ =2, 
postulated by many different particle evolution scenarios, gives 



a = 0.65 all spectra with a < 0.65 should be classified as excep- 
tionally hard. Using this revised criterion, all spectra considered 
here may be regarded as exceptionally hard. 

3.2. Peak levels 

The second relationship between B and 5 derived by Tavecchio 
et al. ( 1998) was obtained from the well-known formula 



syn 



Lc 



(2) 



where U syn and U b are the energy densities of the synchrotron 
radiation field and the magnetic field, respectively, and L s and L c 
are the total luminosities of the synchrotron and the IC emission, 
respectively. The luminosities were calculated from the observed 
spectral slopes and the peak emission levels. Unfortunately, we 
cannot adopt this approach in our estimations because usually 
we do not observe the break in the synchrotron emission (v m ) 
that is related to y m . We instead compare 



MK) = Fc(Vc) 
js(V s ) F s (v s ) 



(3) 



the emissivities (Eq. IA.2I Eq. IB.U with the observed emission 
levels at the peaks. This gives the particle density (Eq. IB.21 > for 
the first part of the particle energy distribution 



(x' s r a 'F c (v c ) 



(4) 



26 where 



. £7**1/3*' 
Jc,sim(x) = — U — I In 



m -n i 



'3x' 

&2 — cc\ \ 4 



( X ') 2 J 

3*' ' 
4y h 



(5) 



is the simplified version of the IC emissivity calculated as a sum 
of two dominant components j CjS im = in n+ in 71 (Tjqs. lBiniB.9t , 
Note that in the component the lower integration bound- 
ary was neglected. Finally, using the synchrotron emissivity (Eq. 
IA.2l i and the transformation given by Eq. IA.5l we can derive an- 
other relation between B and 6 



B = 



\C( ni )(l + z)R 3 K lV ; 
F s (v s )D 2 L 



(6) 



where the luminosity distance (Dl) is calculated in a ACDM uni- 
verse for h = 0.72, Qa = 0.7, and Q,„ = 0.3 for all estimations 
and the spectra calculated in this work. 

3.3. Spectral index in MeV-GeV range 

Another constraint on the basic physical parameters can be ob- 
tained from the slope of the gamma-ray emission in the MeV- 
GeV range. This constraint assumes that the particle minimal 
energy (y m ) is directly related to the radiative cooling inside the 
source. This is one of the main assumptions in our model that we 
discuss in detail in the next section. According to this assump- 
tion 



3m e c 



7m = 



4o- T ct coo \(U B + u m y 



(7) 
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where f coo i is the characteristic cooling time equivalent to the 
source evolution time in this particular case. This and the 
synchrotron energy density (U syD ) are unknown parameters. 
Therefore, for our estimation we assume that U syn = xbIIb is 
energy independent and can be parametrized as some fraction 
(xb) of Ub- In the same way, we parametrize f coo i = xrR/c as the 
multiplied crossing time. Since the spectral peaks are observed at 
almost the same level we may conclude that xb ~ 1. In addition, 
we know that xr must be greater than unity, although, in our es- 
timations we assume that xr > 2. This gives the shortest reason- 
able time for the source evolution. With the above assumptions, 
we can derive an upper limit to the magnetic field value 



B < 



3m e c 2 n 



y 2cr T Ry„ 



(8) 



where y m value required by the above formula can be obtained 
from the spectral slope in the MeV-GeV range. However, there 
is no simple analytic formula that gives y m . 

We derive the the minimum energy (y m ) generating spectra 
from y m = 1 up to y m — yt and searching for the required spec- 
tral index 



a = 



ln{jc(v' 2 )/jc(v[)) 



(9) 



where v\ and v' 2 are the corresponding frequencies to the ener- 
gies 100 MeV and 100 GeV, respectively, in the comoving frame. 
The frequency transformation, which is required for this calcula- 
tion, introduces the Doppler factor dependency into the estima- 
tion. Note that it is necessary to use the full emissivity f Eq. IB. 1 b 
for this particular estimation. The emissivity can be calculated 
for any non-zero value of R, B, and K because the absolute level 
of the emission is unimportant in this calculation. Finally, to cal- 
culate the emissivity we need the energy spectrum break that is 
given by a simple formula 



y b = x c exp 



1 



1 



(10) 



v 2(a 2 - »i) a\ - 1 
derived by Tavecchio et al. 0998). 

3.4. Pair absorption 

The optical depth for the pair absorption inside a spherical 
source can be approximated as 



t(x' ) = -cr T n s (x' )x' 2R, 



(11) 



where n s = 3nj s R/(hcxQ is the number density of soft photons 
per energy interval (e.g. Coppi & Blandford |1990l l. The above 
formula incorporates a simple relation between the soft photon 
and the gamma-ray photon energy x' Q = l/x'. Using the flux 



(12) 



transformation (Eq. IA.5b , we can write 

_ 9 cr T F s (v )D 2 L 
T ~ 10<5 3 (1 +z)Rhc 



which should be smaller than unity for an optically thin source. 
This gives a lower limit to the value of the Doppler factor 



6 > 



9 cr T (l+z) 2a '- l F s (v )D 2 L 



10 



Rhc 



(13) 
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Fig. 2. The upper panel shows the high energy emission of RGB 
J0710+591 and approximated spectra used for our estimations. 
In the lower panel, we illustrate the constraints to the basic pa- 
rameters obtained in four different ways and the values selected 
for the modelling. 



where vo = (m e c /h) / v c . Note that this formula is more accurate 
for a spherical source than the relationship derived by Dondi & 
Ghisellini ( 1995) used in the paper by Tavecchio et al. 0998). 
Our formula infers a lower limit that is almost five times higher 
than the above quoted relationship. 

3.5. Estimation results 

It is difficult to determine the position of the peaks in RGB 
J0710+591 from the available observations. The Swift energy 
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range is quite narrow and the gamma-ray observations made by 
Fermi and VERITAS are not simultaneous. We analysed there- 
fore two alternative solutions focusing first on the Swift and 
Fermi observations and then using the Swift and VERITAS data 
with different assumption about the peak positions. 

To estimate the basic physical parameters of RGB 
J07 10+591 we assume both the position of the synchrotron 
peak (v s = 7 x 10 17 Hz) and the emission level at the peak 
(v s F s (v s ) = 3 X 10~ n erg cm 2 s _1 ). It is also necessary to assume 
the frequency of the IC peak (v c = 9.5 x 10 25 Hz) and the emis- 
sion level at this peak (v c F c (v c ) = 10~ n erg cm 2 s -1 ). However, 
for this peak we assume additional discrepancy of about ±30% 
in both the position and the emission level as well. This discrep- 
ancy is introduced because of the unknown extragalactic absorp- 
tion. The spectral index of the synchrotron emission below the 
peak is assumed to be a\ = 0.5 («i = 2), whereas above the peak 
we use «2 = 3/2 («2 = 4). 

With the above assumptions, we can estimate the values of 
B and 5 in the four different ways described in the previous sub- 
sections. The spectra used for the estimation and the estimation 
result are presented in Fig. [2] On the basis of the constraints ob- 
tained from the peak positions (Eq. []]) and the emission levels 
(Eq. [6jl, we found two crossing areas because of the discrepancy 
in the IC peak. The two other methods provide upper limit to the 
magnetic field value and the lower limit to the Doppler factor. 

To calculate some of the limiting curves, we have to deter- 
mine the radius of the emitting region. This parameter in princi- 
ple can be constrained from the observed variability time scales 
R ct vsa d/(l + Z). However, there is no information about the 
variability of RGB J07 10+591. We therefore selected as large as 
possible a value of the radius R = 2 x 10 16 cm, which still pro- 
vides a good agreement for the four different estimation meth- 
ods. 

In the second approach, we focus on the VERITAS obser- 
vations of RGB J07 10+591 assuming that the IC peak is at the 
frequency v c - 2x 10 26 Hz and the emission level at this peak is 
VcF c (v c ) = 9.5 X 10~ 12 erg cm 2 s _1 . In addition me must assume 
that the synchrotron peak is placed above the observed range 
(v, = 3x 10 18 Hz and v s F s (v s ) = 4 x 10" 11 erg cm 2 s" 1 ). This 
assumption is similar to the solution proposed by Acciari et al. 
( 120101 1 that helps us to explain a "flat" (F ~ 2) intrinsic TeV 
spectrum. 

In the case of 1ES 0502+675, we assumed the same param- 
eters for the synchrotron peak as for RGB J07 10+591 and that 
v c = 2.5 x 10 26 Hz, and v c F c (v c ) = 1 .5 x 10" 10 erg cm 2 s" 1 for the 
IC peak. Moreover, we used ct\ = 0.375 («i = 1.75) and a2 = 2 
(«2 = 5). We explain this particular choice in the next section. 

The main difference between this estimation and the calcu- 
lations performed for RGB J07 10+591 is in the upper limit to 
the magnetic field obtained from the spectral slope in the MeV- 
GeV range. This estimate was calculated under the assumption 
that xb — 1 (t/ S yn — Ub)- However, the observations indicate 
that the emission at the IC peak, in this particular source, might 
even be one order of magnitude higher than the emission at the 
synchrotron peak. This would indicate that xb » 1 but a more 
precise value of this parameter cannot be estimated from the ob- 
servations because the IC scattering at the peak is already in the 
Klein-Nishina regime. This illustrates a limitation of this partic- 
ular estimation method. Finally, we selected as large as possible 
radius (R = 10 16 cm), which provided a good agreement be- 
tween our four different methods of the estimation. 



-9.0 



-9.5 



-13.0 



[ Fermi LAT) 





14 16 18 20 22 24 26 28 
log 10 (i/) [Hz] 




0.15 



CD 

<=* 0.10 



0.05 



0.00 







peaklevels\\ S >12.1, pair limit - 






peak positions 




• (5=13.0, £=0.05 [G] [ : 



10 15 20 
Doppler factor 



25 



30 



Fig. 3. The upper panel the shows multi-frequency emission 
of 1ES 0502+675 observed by the Swift & Fermi instruments 
(Abdo et al. 1201 Obi - squares, Abdo et al. 1201 Pel - circles) and 
the approximate spectra used for the estimation presented in the 
lower panel. The values of B & 6 selected for our modelling are 
indicated by the dot. 



4. Modelling 

The information that we can obtain from the analysed obser- 
vations is significantly limited. We therefore decided to apply 
quite simple model of the emission that is directly related to our 
estimation methods. The basic assumptions of the model that 
we adopted are identical to the assumptions made for the es- 
timations. The VHE originates in a spherical source uniformly 
filled by relativistic particles and an entangled magnetic field. 
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Table 1. The parameters used for the modelling 





n /~< n tat 1 n : rn 1 

RGB J0710+591 


1ES 0502+675 






1st (2nd) estimation 






s 


8(10) 


13 




B 


0.15 (0.13) 


0.05 


G 


R 


2 (1) x 10 16 


10 16 


cm 


Go 


0.25 (3.5) x 10 5 


4.5 x 10" 


cm ' s 


Ji 


0.4 (1) x 10 6 


4.5 x 10 5 




Hi 


3(3) 


4 





The source evolves in time, "fresh" particles are continuously 
injected into the source where they lose energy by means of the 
synchrotron and the IC emission. This is a simple adaptation of 
the jet internal shock model, where the particles efficiently accel- 
erated at a shock front fill the downstream region of the shock. 

The evolution of the particle energy spectrum inside the 
source is described by the kinetic equation 

+ ^ [{f syn + tic } N(y, t)\ = Q(y) + P(y, t), (14) 

where y syn = —Acrjcy 1 UBlOm e c 2 ) is the synchrotron cooling 
rate and the IC cooling rate (yic) is calculated according to the 
prescription of Sauge & Henri (j2004> for an isotropic distribu- 
tion of soft photons, using the full Klein-Nishina cross-section 
in the head-on approximation (Jones 1968). The injection rate is 
described by a power-law function 

Q(7) = Qoy- ni for y>y ; , (15) 

where y, » 1 describes the minimum energy of the injected par- 
ticles. Finally, we simulate the pair creation process within the 
source, where the pair injection rate P(y, t) is calculated using 
the approach described by Sauge & Henri 2004. This process 
in principle should be negligible when we use the parameters 
obtained from our estimations. However, in the estimation our 
intention was to use as large as possible a radius that gives a 
relatively small value of the Doppler factor. The value is only 
slightly larger than the lower limit obtained from our simple es- 
timation (Eq.[T3l). Therefore, the pair creation process may have 
a small impact on the evolution of the particle energy spectrum 
in our more precise simulation. To solve the kinetic equation, we 
used the numerical method described by Chiaberge & Ghisellini 
1999 Calculating spectra for each time step, we are carefully 
checked the energetic balance between the injected and radiated 
energy. 

Such simple time-dependent modelling has been already 
proposed many times to explain evolution of the TeV blazars 
emission (e.g. Dermeret al. 1999, Kusunoseet al. 2000, Bottcher 
& Chiang [20^21 Sauge & Henri [2004b . The main difference be- 
tween our approach and the other models is in the minimum 
injected energy, where in the other approaches y, = 1 is usu- 
ally assumed. In other words, we assume that the acceleration 
process at the shock front is very efficient and "pushes" almost 
all particles towards high relativistic energies. This is somehow 
similar to the formation of either a thermal or quasi-thermal par- 
ticle energy distribution (e.g. Schlickeiser 1984, Katarzynski et 
al. l2006bl) . 

To obtain an exceptionally hard spectrum in the MeV-GeV 
range, we needed to apply a cut-off to the low energy part of 
the spectrum. This cut-off can be produced in a natural way by 
the radiative cooling. Injecting only high energy particles, we 
may obtain a broken power-law particle distribution with the 



break at the y, energy. The minimum energy (y,„) that the par- 
ticles can reach in the source evolution time is given by Eq. [7] 
(e.g. Kardashev |1962| ). This shows that in principle, the particles 
might cool down to y ~ 1, if the source evolution is long enough. 
However, the minimum energy is inversely proportional to time, 
which means that the high energy particles (y ~ 10 5 ) lose en- 
ergy more rapidly than the medium energy particles (y ~ 10 3 ). 
Moreover, in the downstream region of the shock there may be 
another acceleration process, for example a turbulent second- 
order Fermi acceleration. This relatively weak process may not 
produce any high energy activity but it can be strong enough to 
compensate at some point for the radiative cooling and to keep 
the minimum energy significantly above the particle rest energy. 
This is, however, a much more complex approach that requires 
an analysis of possible activity. With no information about time- 
dependent flux variations, we simply assume that the source evo- 
lution time is just 2R/c. 

The spectral index of the particle energy spectrum below the 
break will be constant n\ — 2, whereas above the peak the in- 
dex depends on the injection slope tii — n,- + 1 . This, however, 
is true only for the radiative cooling caused by either the syn- 
chrotron emission or the IC emission in the Thompson regime. 
In the Klein-Nishina regime, the scattering efficiency of the high 
energy particles is significantly lower. This modifies the high en- 
ergy part of the energy distribution (e.g. Moderski et al. 2005). 
It is difficult to describe this effect analytically in our particular 
scenario. We may perform only a simple estimation assuming 
that the energy density of the synchrotron radiation field is given 
by 

f/syn-y J W)dV, (16) 

where v\ is the minimum frequency of the scattered photons 
and v' 2 = min[v' x ,3m e c 2 /(4hy)] introduces the KN restriction in 
a simple form (e.g. Chiaberge & Ghisellini I19991 I. Neglecting 
the lower boundary in the above integral and assuming that 
v' 2 - 3m e c 2 1 '(Ahy) and I s oc (y')~ ffl , we may approximate the IC 
cooling ratio as yic oc -y 2 ^ . Hence, the stationary (N - 0) so- 
lution of the kinetic equation is given by N oc y~ 2 ~ ffl , where the 
slope is significantly lower than the value N oc y~ 2 obtained for 
a constant classical synchrotron or IC cooling. This agrees with 
our more precise numerical tests, which indicate that the radia- 
tive cooling dominated by the IC scattering in the KN regime 
can reduce the particle energy slope by a factor < 1/2. This ap- 
pears to be true for the cooled (below the break) and injected 
(above the break) part of the energy spectrum. Since the emis- 
sion of 1ES 0502+675 seems to be significantly affected by the 
scattering in the KN regime, we assume that ri\ = 1.75 for the 
estimation made in the previous section for this source. We note 
that the energy spectrum index less than two helps in addition to 
explain exceptionally hard MeV-GeV spectra. 

Figure[4]shows the modelling of the RGB J0710+591 emis- 
sion. We present the results obtained for two different sets of 
physical parameters derived from the estimations performed in 
the previous section. Note that the IC spectrum obtained in the 
first case is unable to explain the intrinsic shape of the emission. 
In the second case, the intrinsic spectrum is reproduced well but 
the level of the MeV-GeV emission is significantly lower than 
the observed level. However, the Fermi and the VERITAS ob- 
servations are not simultaneous, hence there may be a difference 
in the emission levels. The detailed values of the physical pa- 
rameters used for the modelling are given in Tab. [TJ There are 
only six important parameters illustrating the simplicity of the 
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Fig. 4. The modelling of the RGB J07 10+591 spectral emission. 
The upper panel shows the first approach where we considered 
the Swift and Fermi observations only. In the lower panel, we 
demonstrate an alternative approach that can explain the ori- 
gin of the TeV emission observed by VERITAS. Note that the 
Fermi and the VERITAS observations were not simultaneous. 
To obtain the observed IC spectra, we used the EBL lower-limit 
model proposed by Kneiske & Dole (2010), where the dashed 
lines show the intrinsic unabsorbed emission. 



model. In addition, we assumed that the maximum Lorentz fac- 
tor of the injected particles is y > 10 7 and the source evolution 
time f evo = 2R/c. 

For 1ES 0502+675, we do not attempt to reproduce the Swift 
UVOT observations. This part of the emission might be domi- 
nated by either the host galaxy or more extended jet structures as 
demonstrated for example in the case of Mrk 501 (Katarzynski et 
al. 1200 11 1. Moreover, no direct correlation between X-ray and op- 
tical variability in other TeV blazars also implies that the optical- 
UV emission has a different origin. Our modelling is also unable 
to explain the inverted spectral slope in the 100 MeV to 1 GeV 
energy range. This is an unexpected and puzzling feature of the 
emission that will require additional investigation. However, this 
may be simply the IC emission produced by the extended jet 
structures, which also dominate the emission in the optical-UV 
range. 
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Fig. 5. The spectra obtained from the modelling of 1ES 
0502+675 for the physical parameters obtained from our esti- 
mation. 



5. Summary 

Several blazars observed by the Fermi Gamma-Ray Space 
Telescope have exceptionally hard spectra in the MeV-GeV 
range. We have demonstrated that in principle all sources with 
a spectral index a < 0.65 should be classified as sources with 
an exceptionally hard MeV-GeV spectrum. Therefore, we have 
attempted to reproduce these spectral slopes by assuming a low 
energy cut-off in the particle energy distribution. This approach 
had been previously developed to explain hard spectra in the TeV 
range and at least in one case (1ES 0229+200 Tavecchio et al. 
120091) it was found to be correct. 

We derived four different methods to estimate the basic pa- 
rameters of the emitting region. Three of these methods had been 
previously proposed and we simply adapted them to the case of 
the low energy cut-off. Our method that gives an upper limit to 
the magnetic field, derived from the spectral slope in the MeV- 
GeV range, was proposed here for the first time. This particular 
constraint can be very useful but it requires accurate measure- 
ments of the spectral index in the not-absorbed MeV-GeV range. 

Finally, we used our modelling estimates that succeed to ex- 
plain the exceptionally hard spectra in terms of particle energy 
evolution. Our model is relatively simple but more complex, and 
hence more accurate modelling would require far tighter obser- 
vational constraints. We have attempted to explain the emis- 
sion of two different TeV blazars. The particle evolution in 
RGB J07 10+591 seems to be dominated by synchrotron cool- 
ing, whereas in 1ES 0502+675 IC cooling appears to dominate. 
In both cases, the estimated parameters of B & 5 provide good 
results in our modelling. This type of modelling is crucial for un- 
derstanding the nature of the blazar emission process, especially 
in non-standard cases such as the exceptionally hard MeV-GeV 
emission. One interpretation of exceptionally hard spectra is that 
for some reason only high energy particles can produce the X- 
ray and gamma-ray activity of TeV blazars. 



Appendix A: Approximation of synchrotron 
emission 

For a power-law particle energy distribution 

N{y) = Ky-\ [cnT 3 ] (A.l) 
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the synchrotron emissivity can be approximated by a simple 
power-law function 

j s (v') = — C(n)KB a+ \VY a , [erg s" 1 cm^sterad^'Hz" 1 ], (A.2) 
4tt 



where a = (n - l)/2 and 
3e 



coo 



4V37re 3 



m e c 2 \ 2nm e c 



r / 3n-l \ r /3//+ \9 



12 



12 



x rfc^|/(8-0fo+l)IX(n + 7)/4)) 



(A.3) 



(e.g. Ginzburg & Syrovatskii 1965 ). Assuming an optically thin 
source, we can write the synchrotron luminosity L s = 4nVj s 
[erg], where V is the source volume. If our source is spherical 
then the intensity of the surface emission is given by 



(A.4) 



where R is the source radius. This gives L s = 4tt 2 R 2 I s and us- 
ing the standard luminosity-to-fiux conversion, we can obtain 
the observed flux F s = n(R/D£) 2 I s [erg cm -2 ], where Di is the 
luminosity distance. For relativistically moving sources at cos- 
mological distances, we have to apply additional transformations 
to the flux 



L 

as well as the frequency v = V 6/(1 + z). 



(A.5) 



Appendix B: Approximation of inverse-Compton 
emission 

The inverse-Compton emissivity can be approximated by a sim- 
ple formula 



jc(x 



N 



3x' 




(B.l) 



where x' = hv' /m e c 2 is the inverse-Compton photon energy nor- 
malized to the electron rest energy and x' is the energy of a syn- 
chrotron photon also divided by m e c 2 . The particle energy spec- 
trum, in our particular estimations, is assumed to be a broken 
power-law 



N(y) = 



Kiy "> for y m <y <y b 
K 2 y~" 2 for y b < y < y x 



(B.2) 



where K 2 - Kiy" 2 ~"' and y m » 1. For such an energy spec- 
trum, we can approximate the synchrotron intensity by a double 
broken power-law function 



f Soix^-oo 
Slix')-" 2 



for x',<x' < x' m 
for x' m < x' Q < x' s 
for x' < x'„ < x' 



(B.3) 



that extends from some minimal energy (jcJ) up to the maximal 
energy (x' x ) and contains breaks at the characteristic energies 
(x' m , x' s ) related to the particle energy by 



r 3.7 x lO^B 



(B.4) 



According to the Eqs lA.2l and lA.41 the normalizing coefficients 
in the synchrotron intensity are given by 



Sl/2 



So 



— C(n 1/2 ) f 

sx(x' m r^ 



RK l/2 B a ^ 



(B.5) 



(B.6) 



Note that inside the source the intensity of the synchrotron radi- 
ation field is on average I s ^ (3/4)/? j s (e.g. Kataoka et al. 1999). 
Therefore, there is no constant 4/3 in the above formulae. 

Since the energy spectrum (Eq. IB. 21 ) and the radiation field 
spectrum (Eq. IB.3l ) are divided into the simple power-law func- 
tions, we can split the inverse-Compton emissivity into six com- 
ponents where each component is described by 



j(a,b)( X 



-\-a b 



(B.l) 



where index a (equal either to 1 or 2) describes parts of the parti- 
cle energy spectrum and index b (equal 0, 1 or 2) indicates parts 
of the synchrotron spectrum. Depending on the a, b values, this 
gives two different results. For a — b,we have 



• t >\ °" r v l 3x ' 

J(a,a)(x ) = — SaK a \ — 



\n(x 2 lxi) 



and for a + b 



j(a,b)(x) = -£-S b K a \^- 



3x'Y a ° x\°-x\» 



a a - a b 



(B.8) 



(B.9) 



Finally, we have to define the integration boundaries {x\,x 2 ), 
which are different for each component. In general, we can de- 
fine the lower integration boundary as a maximum of two ener- 
gies of synchrotron photons 



x\ — max(x' , x' ), 



(B.10) 



where x' p is equal to x' v x' m , or x' s for b - 0, 1,2, respectively, 
and x' q = (4/3)y 2 x' Q is calculated for the characteristic particle 
energies y b , y x for a — 1 and 2, respectively. The upper boundary 
condition is calculated to be the minimum of three energies 



x\ 



mm | x„, x,„ 

p< r 4x , 



(B.ll) 



where x' p is equal to either x' m ,x' s , or x' x for b = 0, 1,2, respec- 
tively and x' q is calculated for y m , y b for a — 1 and 2, respectively. 
The last part of this maximum comes from the simple approxi- 
mation of the Klein-Nishina limit where 



(B.12) 



An example spectrum of the approximate IC emission and a de- 
scription of the integration boundaries and the components of the 
spectrum is presented in Fig. IA.1I This particular spectrum was 
calculated for a set of the parameters that gives all six compo- 
nents of the spectrum. Note that for the parameters that are typi- 
cal for Te V blazars, the IC scattering occurs predominantly in the 
KN regime and only two components (ja,i)) an d 0(2,1)) creates 
the spectrum if 1 << y m « y b . In the case where 1 << y m y b , 
the first part of the energy spectrum (a = 1) and the second part 
of the radiation field (b = 1) basically do not exist. In such a 
case, the spectrum is dominated by the jf^o) component with the 
spectral index in the MeV-GeV range ao - -1/3. However, this 
is an extreme limiting case, never directly observed. 

Acknowledgements. We thank the anonymous referee for a careful review and a 
helpful report. 
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Fig. A.l. An example of an approximation of SSC emission. The figure shows a broken power-law {yt) electron spectrum with a 
low energy cut-off (y m ) that extends up to some maximal energy (y x ). This distribution generates synchrotron emission with two 
breaks (v' m ,v' s ) and extends from some minimal (vj) to some maximal frequency (v' t ). The energy spectrum is divided into two 
parts, whereas the synchrotron radiation is divided into three parts. This provides six different components of the inverse-Compton 
emission. Each component must be integrated with different boundary conditions. The upper right panel demonstrates how to choose 
these conditions. This panel shows synchrotron frequencies as functions of the IC photon frequencies, where the horizontal lines 
indicate the characteristic frequencies in the synchrotron spectrum. Three slanted lines are calculated according to the well-know 
relation vj c = (4/3)y 2 v' :iyn for the three characteristic particle energies {y m ,jb, y x )- The single opposite slanted line shows the Klein- 
Nishina restriction according to the approximation x' = 3/(4xi). The crossing lines are produce six rectangle areas that indicate the 
lower and upper integration boundaries for each component. Each component of the IC emission is indicated by a pair of numbers 
where the first number refers to the part of the electron energy spectrum and the second number shows the part of the radiation field 
used to create this component. 
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